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Time and space integrating folded spectrum techniques utilizing acousto-optic devices 
(AOD) as 1-D input transducers are investigated for a potential application as wideband , 
high resolution , large processing gain spectrum analyzers in the search for extra-terrestrial 
intelligence (SETI) program. The space integrating Fourier transform performed by a lens 
channels the coarse spectral components diffracted from an AOD onto an array of time 
integrating narrowband fine resolution spectrum analyzers The pulsing action of a laser 
diode samples the interferometrically detected output, aliasing the fine resolution compo- 
nents to baseband , as required for the subsequent CCD processing. The raster scan 
mechanism incorporated into the readout of the CCD detector array is used to unfold 
the 2-D transform , reproducing the desired high resolution Fourier transform of the 
input signal 


I. Introduction 

The Search for Extra-Terrestrial Intelligence (SETI) is an 
ambitious program to examine 773 of the nearest solar type 
stars as well as the entire celestial sphere for any evidence of 
intelligent microwave activity, either beamed at us or leaked 
from a civilization’s internal communications (Ref 1). The 
interstellar distances are so large and the resulting 1/R 2 
attenuation is so great that the receiver sensitivity to a signal 
buried in cosmic background noise must be as large as pos- 
sible. Without a priori knowledge as to the transmitted signal 
structure, ideal matched filters cannot be constructed; how- 
ever, a good assumption is that a transmitted beacon would 
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have a narrow band highly coherent structure in a sea of 
broadband white noise. Thus a narrow bandpass filter, with 
resolution as small as 1-100 Hz, could achieve enough pro- 
cessing gain to detect the presence of an intelligent CW signal. 
Since the preferred frequency for interstellar microwave 
communication could conceivably lie anywhere from 0.1 to 10 
or 100 GHz, then from 10 8 to 10 11 frequency measurements 
must be made for each resolvable spot on the celestial sphere 
observed by the radio telescope. If these measurements were 
performed sequentially with a swept local oscillator tuned 
bandpass filter, then as much as a million years of observation 
time could be required A much more practical approach 
would be to build a wideband, high resolution spectrum 
analyzer with a million or more parallel channels of band 
pass filters, thereby cutting down observation time to on the 
prder of a year. Currently two approaches are being pursued, 
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a digital multi-channel spectrum analyzer (MCSA) (Ref. 1) 
with a tree configuration of DFT and FFT subunits, and an 
analog acousto-optic folded spectrum processor as reported 
in this article 

The advanced state of the art of wideband acousto-optic 
(AO) Bragg cells make them ideal candidates for entering 1-D 
data into an optical signal processing system, and the Fourier 
transforming property of a lens makes the implementation of 
real-time spectrum analyzers with 1000 resolvable frequencies 
routine. In this article we explore some possible approaches 
to performing narrow band time integrating (TI) spectrum 
analysis on each of the 10 3 channels of a space integrating 
(SI) AO spectrum analyzer, thereby obtaining a high resolu- 
tion folded spectrum through time and space integration 
(TSI). 


II. Folded Spectrum Processing 

The concept of the Folded Spectrum was introduced by 
Thomas (Ref. 2), who showed that the 2-D Fourier trans- 
form of a raster scanned signal was an orthogonally raster 
scanned version of the long 1-D spectrum of the original 
signal. This is schematically illustrated in Fig. 1. The folded 
spectrum is ideal for an optical signal processing system since 
it utilizes the two available spatial coordinates to perform a 
high resolution spectrum analysis that would be impossible 
in a 1-D system due to the limitations in space bandwidth 
product of 1-D devices. This processing operation can be 
performed in real time by raster scanning a very long input 
signal s(r) onto a 2-D spatial light modulator (SLM), and 
using the 2-D space integrating Fourier transform operation 
performed by a spherical lens. After the complete 2-D raster 
is formed, coherent light is passed through the SLM and the 
resulting amplitude distribution f(x,n ) is Fourier transformed 
by the lens. The power spectrum | F(u c , iy)| 2 is detected by a 
2-D charge coupled device (CCD) detector array at the back 
focal plane. The orthogonally oriented raster readout of the 
CCD turns the 2-D power spectrum back into a high resolu- 
tion 1-D power spectrum l*S'(/)| 2 , at a rate of up to 60 such 
transformations each second. However, the availability of 
highly mature CCD devices with up to 2000 X 2000 pixels 
can not be fully utilized because of the lower resolution and 
dynamic range of currently available 2-D SLMs. Thus an 
architecture for performing the high resolution folded spec- 
trum operation without the need for 2-D SLMs is highly 
desirable. 

Turpin (Ref. 3) and Kellman (Ref. 4) have suggested such 
an approach utilizing acousto-optic devices (AODs) in an 
interferometric time integrating configuration. In this approach 
a 2-D array of local oscillators with coarse and fine frequen- 


cies in x and y is produced by orthogonal AODs to which fast 
and slow chirps are applied respectively. When a frequency 
component of the input signal matches that of a particular 
local oscillator, a baseband heterodyne component is gener- 
ated, and a time integrated peak is created at that location, 
proportional in amplitude to the amplitude of the correspond- 
ing Fourier component. Since this approach involves inter- 
ferometric detection, a bias term will be generated due to the 
local oscillator reference and a signal dependent bias will occur 
uniformly throughout the detector array. These components 
can be separated by placing the desired heterodyne term on a 
spatial carrier. However in a white noise environment, as 
envisioned for the SETI search with a weak narrow band CW 
peak, the signal dependent bias term will produce a large bias 
because all its frequency components are incoherently detected 
at each CCD pixel, thereby using up the available CCD dynamic 
range. Thus the 2-D time integrating technique for folded 
spectrum processing suffers in a multi-signal environment from 
a large signal dependent bias term, and places too large a 
requirement on CCD dynamic range. 

A hybrid approach that utilizes the 1-D space integrating 
Fourier transformation property of a lens in one dimension 
followed by an interferometric time integrating fine resolu- 
tion spectral analysis can combine the best features of the 
2-D SI and 2-D TI approaches to folded spectrum processing. 
This approach utilizes a standard space integrating acousto- 
optic spectrum analyzer (Fig. 3) to perform the coarse spec- 
tral analysis without introducing a signal dependent bias 
term. The number of resolvable coarse frequency bins is 
limited by the time bandwidth product of the AOD and the 
input apodization to approximately T q B = 1000 The doppler 
frequency shift associated with acousto-optic Bragg diffrac- 
tion produces a linear spatial dependence of the temporal 
frequency at the output spectrum that must be converted to 
baseband at all positions for further time integrating pro- 
cessing on a CCD. Bader (Refs. 5, 6) proposed a way to do 
this with a spatially distributed local oscillator (DLO) refer- 
ence wave which is produced with an auxiliary AOD that 
has a repetitive wideband signal applied to its transducer. 
Interferometric detection of the signal and reference Fourier 
spectra produces narrow band baseband signal components 
at each output location. These narrow band signals can be 
used to produce a time integrating fine frequency analysis 
in the orthogonal direction on a 2-D CCD through the use 
of the chirp transform algorithm by premultiplying each nar- 
row band signal with a slow chirp and then convolving with a 
matched chirp using another AOD. This type of system has a 
less severe signal dependent bias build up problem because 
only signals within a single coarse frequency bin contribute 
any bias to that bin. In this article we will present an alterna- 
tive approach to time and space integrating folded spectrum 
processing. 
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This method, which utilizes pulsed laser diodes, was origi- 
nally proposed by Psaltis (Refs. 7, 8) and is illustrated sche- 
matically in Fig. 2. A very long wideband signal is to be 
optically processed utilizing an acousto-optic device as the 
input transducer. However, since the limited time-bandwidth 
product of the AOD will not allow the entire signal to be 
entered into the AOD simultaneously, we will successively 
slide overlapping portions of the long signal into the AOD 
and freeze the acoustic motion by pulsing the illuminating 
laser diode with a very narrow pulse. For each laser diode 
pulse we have a 1-D spatial representation of part of the 
input signal as an optically coherent wavefront, and on succes- 
sive pulses we will transduce successive portions of the input 
signal into the optical system. Thus we have transformed the 
long signal into a 2-D raster of space (. x ), and time or pulse 
number ( n ). The transformation performed by the AOD 
illuminated with a repetitively pulsed laser diode is given by 


N-l 

f(x,n ) = ^ s(t-x/v) xect (x/ A) 8 (t - nT) (1) 

n - 0 


In this expression v is the acoustic velocity of the traveling 
wave AOD, A is the AOD width, and T is the laser diode 
pulse repetition interval. If A > vT, the raster is oversampled, 
which means that successive raster lines are overlapping. We 
perform the desired long 1-D Fourier transform on the input 
signal s(r), with a 2-D transform on its oversampled 2-D 
raster representation f(x,n). Because of the separability of 
the multidimensional Fourier kernel, a 2-D transform can 
be performed as N 1-D transforms along one axis, followed 
by N 1-D transforms along the orthogonal axis This cas- 
caded system can be implemented by performing a spatial 
Fourier transform with a lens for each of the N laser diode 
pulses, producing an intermediate data set F x (u,n ), which is 
also a space-time raster. This is followed by an array of time 
integrating Fourier transform modules that utilize the second 
spatial dimension.^) to perform the desired accumulation of 
partial products. After a succession of N laser diode pulses the 
full folded spectrum is produced in a 2-D format on a CCD 
detector array. The 2-D spectrum can be very simply unfolded 
back into the desired long 1-D spectrum by utilizing the 
raster readout mechanism incorporated in the CCD itself. The 
resulting system utilizes the mature component technologies 
of acousto-optics, charge coupled devices and laser diodes 
to produce a valuable real-time signal processing module for 
high resolution spectral analysis. 


collimated coherent optical plane wave of wavelength X is 
incident upon the AOD at the Bragg angle, 

0 B = -sin -1 A/2 A 0 ~-A/ 0 /2v, 

in order to maximize the optical power coupled into the first 
order diffracted beam. In this expression f Q is the transducer 
center frequency, A Q = v/f 0 is the midband acoustic wave- 
length, and v is the acoustic velocity. A band-limited signal is 
mixed to the AOD center frequency 7^, producing an octave 
bandwidth signal s(t), that is amplified and applied to the 
piezoelectric transducer, thereby launching an acoustic wave 
replica of s(t) into the photo-elastic medium. The acoustic 
wave creates a traveling wave volume index pertubation 
through the photo-elastic effect, which Bragg diffracts some of 
the incident light into the first order diffracted optical replica 
of a single sideband version of the portion of the signal s(t) 
that is currently within the AOD aperture A . In the frequency 
domain this can be viewed as each Fourier component S(f) of 
the signal s(f) within the AOD bandwidth B = 2/ 0 /3, pro- 
ducing a diffracted plane wave of aperture A propogating at an 
angle 6 d ~ Xf/2v with respect to the optic axis, doppler 
shifted by the frequency f due to the acoustic wave motion. 
This diffracted field propogates through a distance F and is 
incident on a Fourier Transforming lens of focal length F, so 
that at a distance F behind the lens the spatial Fourier trans- 
form of the signal within the AOD is formed. At the back 
focal plane the field observed is 


F(x\t) = f s(t-x/v)e- }2irxx 'fi F dxe' OJt 

•Ja 


g/2n At-x/v) 

X e~ l2 ’ nxx ftJ* dx e /u,f 
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s(f) y 2nft 

X A sine [A ( x'/XF - f/v)] df e ,Lot 
[v S(x'v/XF) e ,2nvx ' t ^ V uf ] 


* A sine (Ax' /XF) 


( 2 ) 


III. Space Integrating Spectrum Analyzer 
With a Pulsed Laser Diode 

The operation of a space integrating (SI) acousto-optic 
spectrum analyzer is explained with reference to Fig. 3 A 


In this expression x is the spatial coordinate along the axis of 
acoustic propogation,x' is the spatial coordinate at the Fourier 
plane, and co = 2nc/X is the optical radian frequency. The 
symbol * represents the convolution operation. We have 
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neglected effects of input apodization, acoustic attenuation, 
anisotropic diffraction and band limitation, as well as other 
higher order effects. The exact doppler temporal phase history 
of each Fourier component is reproduced at the back focal 
plane, but the spatial representation of the spectrum is blurred 
by the Fourier transform of the finite aperture function of 
width A thereby limiting the effective frequency resolution 
to 6/ = \/T = v/A, where T a is the acoustic transit time 
across the aperture. This shows that the number of indepen- 
dent analog samples stored within the AOD at any particular 
instant, which is given by the time bandwidth product T a B, 
also gives the effective number of independent frequencies 
that can be determined in the spatial Fourier transform 
plane. The actual frequency resolution is limited by the AOD 
finite aperture sine convolved with the transforms of the 
optical and acoustic apodization functions, in order to pro- 
duce an overall blur spot Q(x f /\F) which yields somewhat 
worse frequency resolution than 1 /T . Normal operation of an 
acousto-optic spectrum analyzer involves the square law 
detection process on a 1-D array of integrating photosensors, 
which incoherently accumulate charge for a time T* in order 
to improve the signal to noise ratio. 

t' 

F(x') = C \F(x',t)\ 2 dt 

Jo 

t' 

= f \S{x'vl-KF)* Q{x'l'KF)\ 1 dt (3) 

Jo 

The output of this type of spectrum analyzer is a power 
spectrum of the signal s(t) with a resolution limited by the 
blur spot Q(u ), with a normalized spatial frequency variable 
u = x f /XF . For a single sinusoid of frequency f input to the 
spectrum analyzer, s(t) - 2a cos (27 r/7), with a single sided 
temporal Fourier spectrum given by S(f) = a 8 (f- /), the time 
integrated output of the spectrum analyzer is given by 

F(u) = T’ \a\ 2 Q 2 (u-f'/v) (4) 

A heterodyne detection technique can be employed by the 
addition of a simple plane wave reference beam. The instanta- 
neous output of such an interferometric spectrum analyzer has 
a linear spatial dependence of the doppler induced temporal 
frequency that is given by 


where the amplitude and phase of the complex spectrum S{f) = 
IS(/)l e/n ^ are reproduced as temporal modulations of each 
resolved frequency bin which is spatially blurred due to the 
apodization. For the case of the single sinusoidal input tone 
at a frequency f\ with amplitude \a\ and phase £2, the hetero- 
dyne detected instantaneous output is 

F(u,t) = \a Q(u - f/v) e ,2lt f * X e ,u>t + r e /a,f | 2 
= \a\ 2 Q 2 (u-f'/v) +\r\ 2 

+ 2\a\r Q(u -f/v) cos(27 r/'f + £2) (6) 

The temporal modulation on the last term thus reproduces the 
input sinusoid in both amplitude and phase, and could be used 
as the input to a fine resolution time integrating processor if 
we first heterodyne it to baseband. A repetitively pulsed light 
source is a simple and effective means of accomplishing this as 
long as the pulse width r is less than the inverse of the band- 
width of s(t), and much greater than the optical period of 
oscillation. When these conditions are satisfied we can consider 
the laser diode pulsed modulation to be approximated by a 
series of delta functions separated by a time T. 

'-(»)■ E «[^] 

n 

= £ S ^-« r ) ( 7 ) 


This series of delta functions will act as sampling impulses on 
the interferometrically detected temporal modulation present 
at the output plane. It is convenient to consider the input 
temporal frequency to be broken up into a coarse and fine 
part by choosing the closest harmonic of the laser diode pulse 
repetition frequency as the coarse term and the offset from 
this frequency as the band-limited fine resolution compo- 
nent. Thus we substitute /' = k/T+8 into Eq. (6), where 
k = mod 1/r /' is an integer and 8 = res < 1/T. For the 
case of repetitively pulsed laser illumination the heterodyne 
output of the AO spectrum analyzer for a single applied 
sinusoidal tone, is given by 


F(u,t) = | [S(uv)e ,27,vut e ,LJt ] * Q(u) + r e>“ f \ 2 F ( u ’ r) “ [kl 2 Q 2 (« - (k/T+ 8)/v) 

= |S( U v)*G(m)I 2 + \r\ 2 +\r\ 2 +2\a\rQ(u-(k/T+6)lv) 

+ 2 r [|S(av)| cos(27TVMf+£2(m>))] *2(u) (5) X cos(27r(£/r + 8)t + O)] p(t) (8) 
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The fine frequency component 5 will not move the blur spot 
Q(u) by a significant amount and can be dropped from within 
Q(u). The coarse frequency component k/T will be sampled 
by the laser diode pulses at times t = nT resulting in an argu- 
ment within the cosine function of 2t mk, with n and k both 
integers, so that this term can be dropped leaving only the 
fine frequency modulation. Thus for the nth laser diode pulse 
at time t = nT the sampled output can be approximated as 

F(u, nT) = \a\ 2 Q 2 (u-k/Tv)+ |r| 2 

+ 2\a\r Q{u - k/Tv) cos(2n8 nT + SI) (9) 

Thus for a single sinusoidal input the interferometric output 
of the space integrating spectrum analyzer illuminated with a 
series of delta functions can be seen to consist of a signal 
dependent bias term located at the position corresponding to 
the coarse frequency component, a uniform reference bias 
term and a sampled interferometrically detected spot located 
at the coarse frequency locus and oscillating in time at the fine 
frequency 6, with phase £2. This effect of the sampling pulse 
train is usually referred to as aliasing and is considered undesir- 
able. It indicates that the sampling rate is too low for the given 
signal bandwidth, and it results in different parts of the signal 
spectrum folding over on top of each other thereby producing 
invalid signal components. However in this case, the coarse 
frequency channeling of the SI spectrum analyzer effectively 
separates in space the various frequency components that 
could alias with each other. It is important to make sure that 
the sampling rate is fast enough to avoid aliasing of the band- 
limited signals within each coarse frequency blur spot Q(u). 
Since the bandwidth of the signal information within each 
blur spot is increased over 1/T a due to the apodization effects, 
we will assume that an appropriate bandwidth would be 2/T . 
To adequately sample these band-limited signal components 
the Nyquist criterion indicates that we should sample at least 
twice as fast as the highest frequency signal component pre- 
sent; thus the minimum sampling rate should be l/J>4/r 
This indicates that the laser diode should be pulsed at least 
4 times as an acoustic signal transits across the laser beam 
aperture illuminating the AOD in order to produce the appro- 
priately oversampled space-time raster 

IV. Time Integrating DFT Based Fine 
Resolution Processor 

A narrow band fine resolution spectrum analysis opera- 
tion requires the coherent accumulation of data samples for 
a period at least as long as the inverse spectral resolution. 
For this reason an SI AO spectrum analyzer is limited in 
resolution to the inverse aperture time of the AOD In order 
to obtain a higher resolution spectrum analysis with AODs a 
time integrating (Tl) approach must be utilized. Coherent 


accumulations can be carried out on an integrating charge 
coupled device (CCD) detector array as long as the ther- 
mally generated detector dark current and optical bias con- 
tributions do not saturate the desired signal components. 
With typical CCD device defect densities, thermal genera- 
tion rates and signal well capacities, it should be possible to 
integrate for a full second, obtaining 1 Hz resolution capa- 
bilities, and with cooled detectors it would be possible to 
integrate even longer For the SETI program these resolution 
capabilities should be more than sufficient since spectral 
coherence of better than 1 Hz would be very difficult to 
obtain after propogating through the interstellar medium 
and being doppler shifted by both transmitter and receiver 
motions. 

The required operation that needs to be performed on the 
output of each coarse frequency bin of the interferometric SI 
spectrum analyzer with a pulsed LD reported in the previous 
section is a narrow band discrete Fourier transform (DFT). 
The DFT is the discrete version of the Fourier transform as 
is appropriate for sampled data streams such as that produced 
by the interferometric spectrum analyzer illuminated by a 
pulsed laser diode. If we examine the third term in Eq. (9) we 
can see that for each input sinusoid we have a sampled fine 
frequency temporal oscillation s(n) that we want to perform 
a DFT on. The DFT is defined for such a sampled input signal 
s(«) as 

N 

S(m) = ^ s(n) e~ l27rnm / N (10) 

n 

This operation is a vector matrix multiplication, where the 
sampled input signal s(n) is the input vector, the DFT kernel is 
the matrix, and S(m) is the output vector. A possible time 
integrating interferometric optical implementation using a 
pulsed laser diode is shown schematically in Fig 4. The light 
from a repetitively pulsed laser diode is split into two by a 
beam splitter. One arm illuminates an acousto-optic modula- 
tor that has a narrowband signal s(f) applied to its transducer 
sampled by the narrow pulse of the laser diode producing a 
temporal modulation s(nT ). In the other arm, a reference 
AOD with aperture T , equal to the LD pulse repetition interval, 
is schlieren imaged onto the output CCD array where it is 
interferometrically combined with the spatially uniform signal 
beam producing a sinusoidal spatial fringe pattern. With each 
laser diode pulse a new column of the DFT matrix is applied 
to the reference AOD transducer, and propogates across the 
aperture to align with the CCD at which time its motion is 
frozen by the LD pulse, thereby imaging that column onto 
the CCD. Thus the reference signal is analogous to a raster 
scanned version of the DFT matrix, except that along the 
columns a continuous representation is used which is later 
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sampled by the spatially discrete CCD pixel structure. Since 
the successive columns of the DFT matrix represent sinusoids 
of successively higher frequencies the reference waveform is a 
stepped frequency chirp, with the constraint that each new 
frequency phase resets to zero phase each T The reference 
signal starts at frequency / and frequency steps by Af each 
time T which is analytically represented as 

r(t) = £ red X cos [2n(f+nAf) (t-NT)] 

(ii) 

This signal propogates m the y direction at a velocity v, and 
aligns with the AOD and CCD apertures at the times of the 
laser diode pulses as given in Eq. (7). The spatial modulation 
diffracted by the AOD on the nth pulse is given by 

r(y, n ) = r(t - y/v) rect [y/vT] 5 (t - nT) 

= rect [ y /VT] e )2n V +n *f)y/v ( ]2 ) 

which is seen to be a finite aperture plane wave with an angular 
spatial frequency that is stepped linearly with the pulse num- 
ber n. The interferometrically detected intensity pattern on 
the nth LD pulse is found by combining weighted signal and 
reference beams. The signal beam propogates at an angle 6 
with respect to the starting frequency of the reference, result- 
ing in a spatial carrier of spatial frequency a = sin 0 /\. 

I(y, n) = | Gs(nT) e j27T ^ a+ ^ y + b r(y , n) | 2 

= \G\ 2 + \b\ 2 +2Gbs(nT) 

X cos [2n (n A fy/v + ay)] (13) 

By time integrating over N successive LD pulses we obtain a 
discrete summation of the interferometrically detected charge 
distributions due to each pulse. 

/O') = /V(|G | 2 + |fc| 2 ) 

N - 1 

+ 2GZ? ^ s(nT) cos[2n(nAfy/v + ay)] 

n = o 

(14) 

This can be recognized as the discrete cosine transform of 
the sequence s(nT), with an analysis bandwidth of Af, on a 
spatial carrier of frequency y = a + (N - l)A//2v. The carrier 
allows a bandpass filter operation to remove the bias terms, 


and the full complex DFT can be reconstructed by measuring 
the phase angle of the carrier for each resolvable frequency 
component. The CCD is composed of an array of pixels that 
must have sufficient resolution to adequately sample the 
spatial fringe structure under the DFT transform. For a pixel 
size Ay, the number of pixels per cycle of the fringe structure 
is I/ 7 Ay, and should be adjusted to 4 to 8 pixels per cycle. 
The total number of CCD pixels required is the product of the 
DFT space bandwidth product with the number of pixels per 
frequency resolution element, which must be at least one 
cycle Thus M = NAfT/yAy > 4NA/T . 

V. TSI Folded Spectrum Based on Tl DFT 

The characteristics of the Tl DFT processor described in 
the previous section are complementary to those of the SI 
AO spectrum analyzer. The Tl DFT processor is a narrow band 
processor with analysis bandwidth Af ( 1/2 T, and high resolu- 
tion of l /NT The SI AO spectrum analyzer is a broadband 
processor with bandwidth £, and resolution ~2/7" a , which may 
also be limited by the CCD pixel size which has a frequency 
width of vAx/\F Since the Tl DFT processor is a 1 -D sys- 
tem, we can spatially multiplex an array of such processors 
in the orthogonal dimension by using a 2-D CCD detector 
array Since each of the Tl DFT processors in the array needs 
to perform the same narrow band spectrum analysis opera- 
tion, and each signal is brought down to baseband by the 
pulsed laser diode, then they all can use the same reference 
AOD. One way that we can configure a 2-D acousto-optic 
folded spectrum processor is as an additive Mach-Zehnder 
interferometer with an SI AO spectrum analyzer to perform 
coarse frequency analysis in the x dimension in one arm of 
the interferometer, and a spatially multiplexed array of Tl 
DFT processors in the other arm in order to perform fine 
frequency analysis in the y dimension on each coarse fre- 
quency bin. The fine frequency analysis bandwidth Af should 
be equal to the coarse frequency resolution per pixel in order 
to produce a full folded spectrum without any gaps. 

The description and analysis of the TSI AO folded spec- 
trum processor is made with reference to Fig 5, which sche- 
matically illustrates the Mach-Zehnder architecture in the 
upper half, and shows cross sections of the fine and coarse 
arms in the bottom half of the figure. The timing of the 
system is controlled by a 10 MHz stable crystal oscillator that 
is digitally divided by 100 to produce a stable 100 kHz pulse 
train used as the laser diode trigger. The laser diode is biased 
just below threshold and pulsed with a 50 nsec pulse every 
10 /isec The anamorphic Gaussian beam profile of the laser 
diode is collimated by the spherical lens and oriented with 
the long axis along the x direction in order to optimally 
illuminate the entire SI AOD aperture. For the initial single 
tone experiments reported in this article the beam splitter 


234 



reflects only about 10% of the light to the SI arm of the 
interferometer, because the processing gain associated with 
the spatial integration performed by the Fourier lens collects 
enough light to make this arm brighter at the output plane, 
so the optical efficiency can be increased by directing more 
light to the TI arm. For an actual SETI signal environment 
of broadband noise, the beam ratio can best be balanced 
through the use of a 50% beam splitter. In the SI arm a cylin- 
drical lens with power in y focuses the light to a narrow slit 
positioned to coincide with the center of the diffracting 
acoustic beam, and incident on the AOD at the Bragg angle. 
The AODs used in this experiment are slow shear mode 
Te0 2 devices with an aperture T q = 70 jusec and a bandwidth 
B - 40 MHz. The signal s(t) applied to the transducer of 
AOD1 produces a diffracted component which is Fourier 
transformed in x by the following spherical lens, which 
simultaneously recollimates in y. Thus each frequency compo- 
nent of s(t) produces a uniform slit of light at the output 
plane, and the position of this slit varies linearly with the 
input coarse frequency. In the TI arm of the interferometer 
a cylindrical lens with power in x is used to focus the light 
down to a narrow slit along y which is incident on AOD2 
at the Bragg angle. Only a 10 Msec portion of the AOD needs 
to be uniformly illuminated with the narrow dimension of 
the LD elliptical beam profile. The reference signal is gener- 
ated by incrementing a 10 bit digital counter with each LD 
pulse, and driving a 10 bit digital-to-analog converter in 
order to make a 1024 level staircase lasting 10.24 msec. This 
signal is then applied to a phase resetable voltage controlled 
oscillator whose phase is reset each LD pulse producing the 
waveform r(t) of Eq. (11), which is amplified and applied to 
AOD2. The diffracted signal from AOD2 is Fourier trans- 
formed by the spherical lens which recollimates in the x 
dimension and allows a bandpass filter to be placed in the 
transform plane. The filtered light is retransformed in y by 
the following cylinder, producing an image of the signal in 
AOD2 with the appropriate scale so that one 10 Msec sinusoid 
fills the entire CCD aperture H in the y dimension, and is 
uniformly spread out in the x dimension. The final beam 
splitter recombines the two beams at an adjustable angle, 
and the CCD time integrates the interferometric product 
between the two waves for a single video frame time of 16 7 
msec. The CCD vertical synchronization signal is used to 
reset the 10 bit digital counter and begin another folded 
spectrum calculation by time integrating the next 1024 
laser diode pulses. 


photogenerated charge profile detected by the CCD on each 
LD pulse is the modulus squared of the field incident during 
that pulse. Thus the time integrated interferometrically gen- 
erated charge distribution recorded on the CCD after a full 
frame of processing is given by 
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X cos [2n (nAfy/v + ay)] (15) 


The phase reset reference wave starts at frequency / and 
steps through N discrete frequencies separated by A/, and the 
angle of incidence in the y direction © can be adjusted to 
make the resulting fringe structure at a desired spatial fre- 
quency 7 = a + (N - 1) Af/2v. The output pattern consists 
of a uniform bias term due to the reference beam, a signal 
dependent bias term with no variation in the y dimension, 
and the last term which is recognized as the 2-D Fourier 
transform of f(x,n) } which is the desired folded spectrum 
of the long 1-D signal s(t) The folded spectrum term is 
riding on a spatial carrier of spatial frequency y in the y 
dimension so that it can be demodulated from the bias terms 
by an appropriate electronic bandpass filter. For the case of 
a single sinusoidal signal of frequency f' - k/T + 6, amplitude 
\a\, and phase f2, the output becomes 


N - 1 

I{u,y) - £ \b e> 2 *V+nAf)y/v 

n = 0 

+ G e l2n(a+ n v >y X a Q(u - k/Tv) e ,2n6nT | 2 
= N\b\ 2 +NG 2 \a\ 2 Q 2 (u-k/Tv) 


The operation of the processor can be described analyti- 
cally by combining the results of the previous two sections 
The optical field incident on the CCD during the «th pulse is 
the sum of the reference wave acting as a linear spatial fre- 
quency modulation in the y dimension, with the spatially 
Fourier transformed signal within the aperture of AOD1. The 


+ 2G\a\ \b\ Q(u - k/Tv) 

X N sine [NT(5 - A fy/H)] 

X cos (2rtyy + O) (16) 
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This pattern is the single tone impulse response of the system 
and is shown schematically in Fig. 6 It consists of a uniform 
bias term, a signal dependent bias ridge at the coarse frequency 
location u = k/Tv , and the folded spectrum impulse of ampli- 
tude \a\ riding on top of the ridge at fine frequency position 
y = bH/Af . The impulse is on a spatial carrier of frequency 
y in the y dimension, with phase angle r2 which will vary on 
successive frames according to its super fine frequency 
A = Te $ 2 /NT^ ’ ^ us P rocess i n g of successive frames can achieve 
even finer resolution. The resolution of Eq. (16) is given by 
the width of the sine function and is limited by the integra- 
tion time to 1/NT; however, the CCD spatial frequency 
response will apodize the DFT and decrease the achievable 
resolution. The fine frequency analysis bandwidth is seen to 
be Af as y is allowed to vary over the CCD height H 

Experimental results from this processor are shown in 
Fig. 7 for a fine frequency analysis bandwidth Af - 4 kHz. 
The processor was operated with a coarse frequency resolu- 
tion per pixel of 80 kHz, but due to limitations of the elec- 
tronics, the fine frequency analysis bandwidth was limited 
to 10 kHz and the data presented in Fig. 7 is therefore only 
a portion of a full folded spectrum. Cross sections along the 
bias ridge are shown in Fig. 7 (a-d) with different fine fre- 
quency inputs. The bias pedestal drops away near the edge of 
the output trace indicating the region of processor analysis 
H = vT. At a fine frequency difference of 120 Hz the construc- 
tive peak is quite large and well away from the edge. At a 
frequency of 240 Hz the destructive peak has moved over by 
about its width indicating that a resolution of 120 Hz or 
better has been achieved. At a frequency of 600 Hz the peak 
has moved over much farther but seems to be broadening. At 
2100 Hz it is clear that the peak is broadening and decreasing 
in size. This appears to be due to the analog nature of the 
reference waveform generation introducing statistical fre- 
quency variations that produce larger phase uncertainties at 
higher frequencies This problem can be alleviated through the 
use of a digital buffer memory and digital to analog converter 
in order to produce the reference waveform However, this 
digital buffer would need to contain 1024 lines of 1024 words 
each that must be read out at a 100 MHz rate, which is cur- 
rently unavailable to us. For this reason an alternative archi- 
tecture with a much simpler reference function was built for 
comparison purposes. 


VI. Time Integrating Chirp Transform 
on TDI CCD 

A common method of performing spectral analysis with 
analog convolvers is the chirp transform algorithm. This 
algorithm results from the observation that the Fourier kernel 


can be decomposed into a convolution and two products 
through the use of the following identity 

-2 nm = (n-m) 2 -n 2 -m 2 (17) 

This expression can be substituted into Eq. (10), which defines 
the desired DFT, and rearranged to yield the chirp transform 
expression. 
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S(m) = ^ s(h) exp (-/27T nm/N) 

n= o 

iV-l 

= exp (- jirm 2 /N ) ^ [$(/?) 

n = 0 

X exp (-JTrn 2 /N)] exp (Jn (n - m) 2 /N) (18) 

An examination of this representation indicates that we can 
accomplish a DFT by first pre-multiplying the input signal 
with a chirp, convolving the resulting product with a matched 
chirp, and finally post-multiplying the convolver output with 
a matched chirp in the Fourier domain. Often only the power 
spectrum is required in which case we can neglect the final 
post-multiplication by the quadratic phase factor, since its 
magnitude is unity everywhere. 

There are a large number of methods for performing con- 
volutions with time integrating optical techniques utilizing 
acousto-optics and charge coupled devices The features which 
must be incorporated in a convolver are an array of multipli- 
cations, a shift operation and an integration. An architecture 
we have used with great success for chirp correlations (Ref. 10) 
requires a temporally modulated input, a fixed chirp reference 
mask, and a CCD operated in the time delay and integrate 
(TDI) mode, where the detector time integrates in a sliding 
coordinate frame. The CCD readout circuitry is modified so 
that it operates in a continuously scrolling fashion referred 
to as TDI or shift and add In this manner the CCD accumu- 
lates photogenerated charge at each pixel site for one or more 
laser pulses, the charge is then transferred vertically (y) by one 
pixel of width Ay, and each pixel accumulates more charge, 
adding it to that previously detected. By the time a pixel has 
fully traversed the array of height Y, it has accumulated photo- 
generated charge from each position in space at successively 
delayed times If a temporally modulated sampled signal 
f(nT) is uniformly incident on the chirp mask 

c(y ) = cos (by 2 ) rect ( y/Y ) 
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and then imaged on the TDI CCD the resulting signal that 
scrolls off the edge of the detector is given by the convolution 
of the temporal signal with the spatially sampled image of the 
mask 

I(mT) = ^ f(nT) c(mAy) *5 (m-n) 
n 


= f(nT)c((m-n)Ay) (19) 

n 

The index m enumerates the output samples obtained at the 
edge of the CCD at successive times mT. We can use this con- 
volver to configure a time integrating interferometric chirp 
transform power spectrum analyzer by introducing a quad- 
ratically phase modulated reference beam, and interfering 
this with a baseband phase modulated signal beam in order 
to perform an interferometric chirp pre-multiplication of the 
input signal as the detected intensity 

f(nT) = s(nT) cos (h (nT) 2 V 2 ) 

where V = A y/T is the TDI CCD shifting velocity The result- 
ing time integrated power spectrum output would be obtained 
by demodulating the interferometric product term from the 
bias terms and squaring the chirp transform electronically. The 
analysis bandwidth of this type of TDI chirp transform spec- 
trum analyzer is limited by the sampled mask resolution to 
half the CCD shift rate for each sideband of the reference 
chirp. A two sided chirp reference can cover an analysis 
bandwidth equal to the CCD shift rate A feature of this type 
of chirp transform is that the spectral amplitude decreases 
linearly away from the center frequency as does the spectral 
resolution, due to a decrease in the matching chirp overlap. 

The chirp transform spectrum analyzer can be readily 
incorporated in the TSI folded spectrum processor by includ- 
ing a quadratically phase modulated reference beam, and a 
chirp reference mask imaged onto a TDI CCD The resulting 
system architecture is shown schematically in Fig. 8 and can 
be seen to be just a modification of the previous optical 
system. The analytical description of the TSI folded spec- 
trum system with the TDI chirp transform fine resolution 
processors results in the same signal and bias terms as in the 
previous system shown in Fig. 6. However the TDI approach 
has the advantageous feature that it averages coherent optical 
artifacts and CCD pixel noise along the shifting column 
thereby minimizing the system fixed pattern noise. The 
system output is given by 


TV-1 
n = Q 

+ r exp [-jb (nT) 2 V 2 ] 

X cos (b(mAy) 2 ) rect [mAy/Y] *5 (m - n) 

= V* f f f(x,n ) e~ ,2nux dx 

h Ua 

X 2 r cos [b(nAy) 2 ] cos ( b((m - n) Ay) 2 ) 

X rect [((m - n) Ay)/Y] + bias terms (20) 

We can see that the quadratically phase modulated reference 
beam results in an interferometric premultiplication of each 
sampled baseband coarse frequency bin with a temporal 
chirp. The combination of the chirp mask and the TDI CCD 
results in a chirp convolution for each coarse frequency 
channel, and the resulting folded spectrum raster scans out of 
the CCD as it scrolls. Experimental results from this processor 
are shown in Figs 9 and 10. Figure 9 (a-d) shows cross sections 
through the coarse frequency ridge of the impulse response 
of Fig. 6 for several input fine frequencies, and the decrease 
of peak height with increasing frequency is apparent. Figure 10 
shows the 2-D Folded Spectrum output displayed on a video 
monitor for a variety of input frequencies, where the coarse 
ridge is seen to move with 10 MHz frequency change, and 
3 kHz moves the fringes in the fine dimension along the 
ridge The spatial and temporal chirp references used in these 
experiments were two sided with 9 fringes per side giving a 
compression of 36, and resulting in a spectrum analysis with 
36 independently resolvable frequencies. The sideband struc- 
ture is highly visible due to the hard clipping of the mask, 
both the spectral resolution and the linearity are better than 
100 Hz, and the performance closely matches the expected 
theoretical behavior. 


VII. Demodulation of Signal From 
Bias Terms 

A crucial requirement for the successful operation of a 
TSI signal processing system is the ability to demodulate the 
signal from the bias terms. There are two main approaches to 
signal demodulation which are shown in Fig. 11, spatial carrier 
demodulation and CCD based demodulation. In carrier demod- 
ulation the interferometric signal term is placed on a spatial 
carrier by the adjustment of the reference beam incidence 
angle, while the bias terms remain at baseband When these 
signals are output from the CCD, a bandpass filter can sepa- 
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rate the different components, and allow measurement of 
spectral amplitude and phase. The filtered output can be 
electronically squared and low passed in order to obtain a 
power spectrum, or rectification followed by a low pass can 
be used for measuring spectral amplitude components. The 
CCD raster readout mechanism can be used to change a spa- 
tial carrier in the horizontal direction into a fast temporal 
carrier that can be simply bandpassed by an active filter. A 
spatial carrier in the vertical dimension requires a more com- 
plicated set of video delay lines and a tapped filter in order 
to implement the desired bandpass operation. The spatial 
carrier can be placed in the coarse or fine frequency direc- 
tions, but will cut down the available CCD space bandwidth 
product in whichever dimension it is used. In order to avoid 
losing detector resolution, another possible approach is to have 
two synchronously operating CCD detectors, one to detect 
the interferometric signal and bias terms, and the other with 
just the signal term in order to detect the signal dependent 
bias. The two detectors must be precisely aligned and sub- 
tracted pixel by pixel in order to remove the signal dependent 
bias; the reference beam bias term should be uniform across 
the array and can be subtracted with a level shifter. An experi- 
mental bias subtraction technique we are investigating requires 
a special purpose CCD with a fill and spill structure abutted 
to the transfer register in order to bleed off an appropriate 
amount of charge to leave the signal term and a small amount 
of bias, rather than have the bias use up the available CCD 
dynamic range. 

VIII. AO Folded Spectrum Processor 
Performance 

A summary of the performance parameters for the TSI 
folded spectrum processor described in this article are given 
in Table ] , for the Caltech prototype, and for the potential of 
a fully engineered state of the art brassboard. The processor 
analysis bandwidth is determined by the AOD bandwidth, 
the CCD width and Fourier lens focal length, and by the 
laser diode pulse width at an acceptable level of illumination 
and temporal coherence. The Caltech prototype was limited 
by the CCD width, which can easily be modified The mini- 
mum laser diode pulse width is determined by the multimode 
transient operation for the first 1 nsec which results in inco- 
herent operation inappropriate for this processor (Ref. 11); 
however, light limitations will also dictate a minimum pulse 
width determined by peak output power. The number of 
coarse frequency bins is determined primarily by the number 
of CCD pixels, but another limitation is due to the AOD time 
bandwidth product and apodization uniformity. The avail- 
ability of the Tektronix 2000 X 2000 CCD is very encouraging 
for the development of optical processors requiring high 
resolution 2-D detectors. The coarse frequency resolution is 
limited by the AOD transit time and apodization, but also 


determined by the FT lens focal length and pixel width. The 
number of fine frequency resolution bins is determined by 
the space bandwidth of the reference functions and the num- 
ber of CCD pixels available. If a spatial carrier is included in 
the fine frequency dimension, then the number of fine resolu- 
tion bins is further limited, and a maximum spatial carrier 
frequency of 4-6 pixels per cycle indicates that a maximum of 
500 fine frequency bins could be predicted The fine resolu- 
tion is fundamentally limited by the CCD integration time 
during which the LD is being pulsed, and the prototype 
processor has closely approached this limit. Actually a greater 
difficulty is encountered in attempting to achieve a fine 
frequency analysis bandwidth equal to the coarse frequency 
resolution, which requires high bandwidth electronic refer- 
ence generation currently not obtainable with our prototype 
circuitry. The TDI chirp transform approach is limited in 
analysis bandwidth by the mask resolution and by the obtain- 
able CCD shift rate, which is currently limited to 15 kHz for 
video compatible devices, but could be increased with special 
CCD circuitry. The total number of independent spectral 
measurements made in our initial experiments was only 
15,000 but a potential for on the order of 10 6 clearly exists. 
The dynamic range is currently quite small due to coherence 
limitations and light starvation, but this will be greatly im- 
proved in the next generation processor. The most important 
performance parameter for the SETI search is the sensitivity 
to a signal buried in noise, and the processing gain limitations 
of the analog technique used to perform the spectral analysis 
The large AOD dynamic range allows very weak signals to be 
effectively transduced into the optical processor, but a huge 
CCD dynamic range is not required to achieve the desired 
processing gain, since the spatial integration achieves much 
of it. 

Since sensitivity and processing gain are so important to 
the SETI application we will include a simple analysis of the 
processor limitations based on the model presented in Fig. 12. 
For a weak narrowband CW input buried in white noise with 
power o 2 we will have a system input 

5(r) = + n( 0 (21) 

The amplifier and AOD have a large linear region and a satura- 
tion regime that will introduce third order intermodulation 
products that will degrade system performance An ideal gain 
G is included due to an optical power increase which will 
not increase the relative power in the intermodulation terms. 
The processor consists of an SI coarse frequency channelizer 
which acts as 10 3 bandpass filters (BPF), cascaded onto an 
array of 10 3 TI spectrum analyzers, resulting in the synthesis 
of 10 6 narrow band filters. The detected spectral amplitudes 
are corrupted by the inclusion of signal dependent noise, and 
the detector noise term with variance . The processing gain 
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is defined as the ratio of input signal to noise ratio (SNR) to 
output SNR, and ideally it would be given by 


SNR 


PG = 


out 


SNR 


a 2 IP 2 a 2 

« 2 lot 


= 1 IP 2 


10 6 


( 22 ) 


where 0 2 is the BPF width divided by the entire bandwidth, 
which is equivalently the inverse of the number of channels 
After the demodulation operation, which will remove the 
bias terms, including the bias shot noise contribution in 
a 2 , the expected value of the power spectrum component 
at the frequency of the input tone will be 


(S 2 J = a 2 G 2 + /3 2 G 2 o 2 + o 2 (23) 


The input sensitivity is determined by the minimum detectable 
signal at the desired false alarm probability, which for SETI 
is on the order of 10 -10 ~ e" 20 , in order to keep the total 
number of false alarms manageable By cross scanning the 
sky and comparing positions and frequencies of spectral 
components exceeding the given threshold, the false alarm 
probability can easily be squared. Thus the minimum detect- 
able signal resulting in a false alarm probability of e" 10 would 
be 10 standard deviations above the noise 


a 2 G 2 > 10 [P 2 G 2 a 2 + a 2 ] 

mm ir n d 1 


(24) 


To minimize effect of the detector noise limitations we 
increase the gain G to the point that the signal dependent 
noise and the detector noise are comparable. In this case 
minimum input SNR that can result in a signal exceeding the 
threshold is given by 


SNR . 
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a 2 . /o 
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= 10 [J3 2 + a 2 /G 2 a 2 n ] ~ 20/3 2 (25) 


Thus the resulting input sensitivity is degraded by 13 dB from 
the initial -60 dB, but only 3 dB of the degradation are an 
artifact of the analog optical processor. This estimate for the 
actual input sensitivity of -47 dB will be further degraded by 
the intermodulation product terms due to the amplifier and 
AOD nonlinearities, by the shot noise due to the bias, and 
probably most importantly by light starvation limiting the 
gain G 


IX. Conclusion 

In this article we have presented preliminary experimental 
results from time and space integrating folded spectrum pro- 
cessors that may have potential application in the SETI 
program. The sampling action of the repetitively pulsed laser 
diode eliminates the need for a distributed local oscillator 
reference wave by aliasing all of the fine frequency compo- 
nents to baseband at the output of the coarse frequency 
channelizing space integrating spectrum analyzer as is required 
for the following CCD based narrow band time integrating 
spectrum analyzers. Two methods of time integrating fine 
frequency analysis were presented, the time integrating DFT, 
and the time delay and integrate chirp transform. The TI DFT 
has the attractive possibility for incoherent addition of spectra 
in order to improve SNR, but requires a very complex and 
accurate reference waveform. The TDI chirp transform helps 
average out coherent noise and detector noise, and has a 
simple reference function, but does not allow simple complex 
spectral demodulation, and is bandwidth limited by the 
available CCD shift rate The possible advantage of this tech- 
nique over alternative digital approaches is a decrease in sys- 
tem cost and complexity, along with the elimination of the 
high speed analog to digital converter required for digital 
processing. For a possible future spaceborne SETI search the 
optical processor has decreased size, weight and power require- 
ments over that of a digital approach 
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Table 1. AO processor performance 



Parameter 

Caltech 

Prototype 

Potential 

B = 

Bandwidth 

30 MHz 

1 GHz 

"c = 

Course Bins 

384 

2000 

II 

Course Resolution 

80 kHz 

5/1000 

7 = 

Fine Bins 

36 (actual) 
491 la 

carrier 

(potential) 

2000/a 

' carrier 

II 

Fine Resolution 

120 Hz (actual) 

1 kHz (potential) 

fj 1000 

jV = 

Spectral Resolution 

15,000 

10 6 

DR = 

Dynamic Range 

50 levels 

1000 levels 


Sensitivity 

-40 dB 

-60 dB 
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RASTER 




Fig. 1. Two-dimensional space integrating folded spectrum pro- 
cessing: (a) the relationship between a long 1-D signal, its high 
resolution spectrum, the folded spectrum, and the signal’s raster 
representation; (b) a space integrating folded spectrum processor. 



F(u, v) = Yl [I n ) ex P H2ttux) dx] exp (-j27rvn/N) 
n=0 J 


Fig. 2. Time and space integrating folded spectrum processor showing the conversion 
to a space-time raster, spatial Fourier transformation cascaded onto an array of 
temporal Fourier transforms. 
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Fig. 6. Impulse response for TSI folded spectrum processor for a 
single tone input at frequency f ’ = k/T + 8, showing uniform bias, 
signal dependent bias, and impulse riding on carrier 


(c) Sf = 600 Hz 


Fig. 7. Experimental results for the TSI folded spectrum 
processor of Fig. 5, showing cross sections through the 
coarse frequency ridge for various fine frequencies 



244 




SI REFERENCE CHIRP CHIRP TDI 
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Fig. 8. Architecture for performing TSI folded spectrum using a TDI CCD to perform a 
fine resolution chirp transform 
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(a) 8f = -100 Hz 


(b) 8f = -1000 Hz 


(c) 8f = -2000 Hz 


(d) 8f = -3000 Hz 



Fig. 9. Cross sections through coarse frequency ridge for the folded 
spectrum processor of Fig. 8 for various fine frequencies 


246 



«ero0NAT_M|os 

OF PO0R~ D l,r: "TY 

ORIGINAL PAGE IS 
OF POOR QUALITY 


(a) 60 000 MHz (b) 60 003 MHz 




(c) 50 000 MHz (d) 50 003 MHz 

Fig. 10. Impulse response of TSI folded spectrum processor displayed in a two-dimensional format on a video monitor 
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(a) CARRIER DEMODULATION 

SIGNAL BANDPASS RECTIFICATION ENVELOPE 

+ BIAS FILTER DETECTION 





SIGNAL 

DEPENDANT 



Fig. 12. Model used for sensitivity and processing gain analysis 


(d) FILL AND SPILL ADAPTIVE BIAS SUBTRACTION CCD 



Fig. 11. Demodulation schemes required to separate signal 
from bias terms 
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